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SUMMARY 


In this paper, the Boundary Force Method (BFM) , a form of an indirect 

boundary element method, is used to analyze composite laminates with cracks. 

The BFM uses the orthotropic elasticity solution for a concentrated horizontal 

and vertical force and a moment applied at a point in a cracked, infinite sheet 

as the fundamental solution. The necessary stress functions for this 

fundamental solution were formulated using the complex variable theory of 

orthotropic elasticity. The current method is an improvement over a previous 

method that used only forces and no moment. The improved method was verified by 

comparing it to accepted solutions for a finite -width , center-crack specimen 

subjected to uniaxial tension. Four graphite/epoxy laminates were used: 

[0/+45/90] , [0], [±45] , and [±30] . The BFM results agreed well with accepted 
s s s 

solutions. Convergence studies showed that with the addition of the moment in 
the fundamental solution, the number of boundary elements required for a 
converged solution was significantly reduced. Parametric studies were done for 
two configurations for which no orthotropic solutions are currently available: a 
single edge crack and an inclined single edge crack. 
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INTRODUCTION 


Previous work for orthotropic materials (ref. 1) used the stress functions 
for a concentrated horizontal and vertical force applied at a point in a 
cracked, infinite sheet as the fundamental solution. The present work extends 
the fundamental solution to include a moment, as in the method for isotropic 
materials. The necessary stress functions for this fundamental solution are 
formulated using the complex variable theory of orthotropic elasticity (ref. 2). 
The addition of the moment DOF to the orthotropic formulation of the BFM is 
evaluated by comparing the BFM results to accepted solutions for a finite-width 
center-crack specimen subjected to uniaxial tension. Convergence studies are 
done comparing the methods for forces only and for moments and forces. 

Parametric studies are done for specimens with a single edge crack and an 
inclined single edge crack. 
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half length of a center crack or length of an edge crack, m 
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constants in stress functions, N/in 
complex constants (i,j = 1,2) 

Young's moduli in the x- and y-directions , respectively, Pa 
mode I and mode II stress - intens i ty correction factors 
influence coefficient matrix, N/m 
orthotropic shear modulus, Pa 
height of plate, m 

mode I and mode II stress - intensity factors, Pa,/m 
unit moment on the i^ boundary element, N-m 
concentrated moment, N-m 
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number of boundary elements 

normal and shear crack- face loading functions, Pa 

units loads on the i^ boundary element, N 

concentrated forces in the x- and y-di rections , 
respectively, N 

vector of unknown forces, N 

x- and y-components of applied loading, N 

vector of external loads, N 

remote applied stress, Pa 

location of load point on crack face, m 

externally applied moment, N-m 

complex variable, (i = 1 to 4) 

width or half -width of plate, m 

Cartesian coordinates, m 

load point, (z Q = x Q + iy Q ) , m 

incremental distance, (i = 1 to N), m 

orthotropic stress functions, Pa 

Poisson's ratios 

stresses, Pa 

roots of the characteristic equation 


BOUNDARY FORCE METHOD 


The Boundary Force Method (ref. 3) is a numerical technique which uses 
fundamental solutions for concentrated forces and moments in an infinite sheet 
to obtain the solution to the boundary value problem of interest. These 
fundamental solutions are used to exactly satisfy the stress-free conditions on 
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the crack faces. The other boundary conditions are approximately satisfied by 
applying the appropriate sets of concentrated horizontal and vertical forces and 
moments along the boundary. 


Consider, for example, the finite cracked plate subjected to uniaxial 
tension shown in Figure 1. In the BFM, an imaginary boundary corresponding to 
the finite Pi ate is traced on a cracked infinite sheet. These boundaries are 
then divided into a finite number of boundary elements. On each boundary 
element, a concentrated force pair and and a moment fh (i — 1 to N) 

are applied at a small distance S ^ on the outward normal from the mid-point of 
the boundary element. This small offset from the boundary element was used to 
avoid the inconvenience of singularities in the computation of the stresses on 
the boundaries. In the present work, 6^ was set equal to one quarter of the 
element length. 

til 

The concentrated forces and and moment acting on the l 

boundary element produce stresses on the other boundary elements. The resultant 
forces and moments on each boundary element j are found by integrating the 
stresses over the boundary element length (assuming a unit thickness) . On each 
boundary element j , the sum of all the resultant forces and moments must be 
equal to the externally applied forces and moments acting on the boundary 
element. Denoting the externally applied forces as R and R , and the 

X j y j 

externally applied moment as T^ , we can write the following equilibrium 
equations for the j ^ boundary element: 



(F P. + F Q. + F M.) 
, x.p. l x.q. l x.m. l 
1 J 1 J 1 J i 
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( 1 ) 


R - E (F P.+F Q. + F M. ) 
y. . , y.p. i y.q. 1 y.m. l 

y 2 1=1 -^J^i ^J H i i 


T. = £ (C. P. + C. Q. + C. M.) 

J i= i JPi 1 1 J m i 1 


Here N is the total number of boundary elements; F , F , F , F , 

J x.p. y-p. x.q. y.q. 

J i i J i J J M i 

F . F , C. , C. , and C. are called influence coefficients and are 
x.m. y.m. ip. iq. im. 

j l J j l Jt i JM i i 

defined as follows. 


F - force per unit force in the x-direction on the j boundary 

X j P i 


element due to unit load p^ acting in the y-direction on the i 
boundary element. 


. th 


F = force per unit force in the y-direction on the j t ^ 1 boundary 

y j P i 


element due to unit load p^ acting in the y-direction on the i 
boundary element. 


.th 


fcl"l 

C. = resultant couple per unit force created on the j boundary 

JPi 


element due to unit load p^ acting in the y-direction on the i 
boundary element. 


. th 


, F 


x . ci . 1 y.q.’ x.m.’ y.m.’ j q . 

j * l j l l i 


C. , and C. are defined in a similar manner. 


j™ 
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The resulting system of equations for N boundary elements can be written 


as 


^3Nx3N {P, 3Nxl “ (R, 3Nxl (2) 

where [F] is the influence coefficient matrix, { P } is the vector of 
unknowns, P, Q, and M, and (R) is the vector of externally applied forces. 
The influence coefficient matrix is square, fully populated and non- symmetric . 

Because the influence coefficient matrix and the externally applied load 
vector are known, the unknown force vector can be obtained by solving the system 
of linear algebraic simultaneous equations. The calculated set of P Q. , and 
th acting on the imaginary boundaries in the sheet will approximately satisfy 
the required boundary conditions, and, thus, produce a stress distribution 
inside the imaginary boundaries that is approximately equal to the stress 
distribution of the desired boundary value problem. 

DERIVATION OF FUNDAMENTAL SOLUTION FOR ORTHOTROPIC MATERIALS 

For orthotropic materials, the BFM uses the elasticity solution for a 
horizontal and vertical concentrated force and a concentrated moment in an 
infinite orthotropic sheet with a crack. The formulation of this solution is 
presented below. 


Stresses in Orthotropic Materials 

From Lekhnitskii (ref. 2), the stresses in an infinite orthotropic sheet 

t t 

can be written in terms of the two stress functions, <j>^ and <f>^, as follows: 
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a x ( x > y) = 2 Re 


2 ' 2 ' 

(z i > + ^ 2 ^2 ^ z 2^ 


® (x,y) = 2 Re 


( Z l> + ^2.^ Z 2 ) 


( 3 ) 


xy 


(x,y) 


-2 Re 


M 1 ^ 1 (z 1 ) + n 2 <f> 2 (z 2 ) 


For simplicity, the derivation of the stress functions for the concentrated 
forces and for the moment acting at an arbitrary point in a infinite, cracked, 
orthotropic sheet will be presented separately. 

Concentrated Forces . The derivation of the stress functions for a pair of 
concentrated forces acting on a cracked, infinite orthotropic sheet was 
presented in ref. 1. For completeness, the derivation of these stress functions 
is shown in Appendix A. The stress functions are given below. 




r (A * 2 * 
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z, - w, f ( z l> w i> a ) + T f( z i> w ]_> a ) 
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Z 1 * *2 


Z 1 ' W 1 
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7 2 ( z 2 ) 


2(Mi - An) J 
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^ ' ■» 2 
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(Mi - 
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where f(z, w, a) - J z 2 - a 2 - ./w 2 - a 2 - z + w, and the other terms in 
equations (4) are defined in Appendix A. 


Moment . Figure 2 shows how superposition was used to determine the 
stress functions for a moment acting on a cracked, infinite orthotropic sheet. 
The uncracked sheet with the moment applied at z^ is shown in Figure 2(b). 
The normal and shear stresses acting on the line y = 0, |x| < a, are shown as 

N(x) and T(x) . The problem with the crack-face loading (Figure 2(c)) is 
superimposed on the uncracked problem (Figure 2(b)) to produce the stress-free 
crack face shown in Figure 2(a). 


The stress functions for a moment in an infinite, uncracked orthotropic 
sheet (used for Figure 2(b)) are derived (following ref. 4) in Appendix B. 
These stress functions are 


-A 


M 


^l (z l ) ” (Zj_ - w L ) 2 


. ' , . ’ B M 

W ’ (z 2 - w 2 ) 2 


(5) 


where 


1 C 1 1 ' C 12^ M 


b m = 


(a* 2 C 2 1 * C 22^ M 


The terms , C^, C^. C 2 \’ anci C 22 are cle ^^ ned *- n Appendix A. 


The stress functions for the loads N(x) and T(x) applied to the crack 
face in Figure 2(c) are given in Appendix A, equations (A4) . Here, 
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N(x) - a (x,0) and T(x) = r (x,0) in the uncracked sheet in Figure 2(b). 

y *y 

From equations (3) and (3), 


ff (x, 0) = -2 Re 


r (x,0) - 2 Re 


r + — fM — l 

L(z 1 - w 1 ) 2 (z 2 - w 2 ) 2 J 


f M 1 A M M 2 B M 1 

L(z 1 - w 1 ) 2 (z 2 - w 2 ) 2 J 


(6) 


Using equations (6) and remembering that 2 Re[f(z)] = f(z) + f(z) , the 
loading functions in the integrands in equations (A4) can be simplified as 
follows : 


VSl 


P 2 B M 


P 2 N(x) + T(x) - -2 ^Re [ (x -"'ji + (x 1 + 2 Re [ (x - w^ 2 + (x - « 2 ) 

^2 ' "l )A M (tl 2 ' K )k H ( ^2 ~ ^2 )B M 


(x - w ]L )‘ 


(x - w 1 ) i 


(X - w 2 )‘ 


(7) 


^N(x) + T < x > “ - 2 ^l Re t(x - Wl ) 2 + (x - w 2 ) 


B M , , M 1 A M M 2 B M 

2 +2 Re [ (x _ )2 + (x . )2 . 


( ^1 - V A M _ °*1 ' ^2 )B M _ (/J l 1 m 2 )B M 

(x - w x ) 2 (x - w 2 ) 2 (x - w 2 ) 2 


Substituting these expressions in equations (A4) , integrating, and then 
adding equations (5), the following expressions are found for the stress 
functions for the loading shown in Figure 2(a): 
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-1 


where 




2(^1 - Hy) Jz \ - 


[ (m 2 - M 1 ) A m D(z 1 ,w 1> a) 


^ ( M 2 " M j ) A^j U(z 1 ,w L ,d) 4 (^2 


/‘ 2 ) D ( z j . w 2 .a) 


-A 


M 


(Zl - Wj ^) 2 


^ 2 (z 2 ) 


2 (Mi - Mo) / 


ZA - a‘ 


[ (M x - M 2 ) B m D(z 2 ,w 2 ,a) 


1 v 2 
+ (m x - M 1 ) A m D(z 2> w 1 ,a) + (m-l 

. ~ b m 


M 2 ) B m D(z 2> w 2 ,a) ] 


(z 2 • W 2 )2 


( 8 ) 


_ . . (w^ - a^ 1 ,/z 2 - a 2 + (a^ - wz) ,/w 2 - a 2 

D(z,w,a) - ~ ~ o 

(w - a )(z - w) 


The stress functions used in the fundamental solution for orthotropic 
materials are found by combining equations (4) and (8). Once the stress 
functions are known, it is a simple step to calculate the stresses at any 
point in the body using equations (3). Then, as mentioned earlier, the 
stresses are integrated over each element length to obtain resultant forces 
for use in equation (2) . 

Stress - Intensity Factor Equation 

From Snyder and Cruse (ref. 5), the stress - intensity factors for 
orthotropic materials may be expressed in terms of the stress functions as 
follows : 
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= -2j2n 


(9) 


K T + 

I M 


II 


**1 ~ "2 


lim { J z, 


4 > (z^) ) 


By substitut ing from equations (4) and (8) into the above equation and 
taking the limit, the mode I and mode II components of the stress- intensity 
factor for a horizontal and vertical force and a moment in an infinite cracked 
orthotropic sheet may be written as follows: 


K 


I 




+ ( A * 2 


( A G < w 1 ,a) + ^ H(w 1 ,a)] 
/i 1 ) [A G(w 1 ,a) + A^ H(w^ , a) ] 


+ (/i 2 - p 2 ) [B G(w 2 ,a) + B m H(w 2 ,a)] 


( 10 ) 


where 


G(w, a) 


./w 2 - a 2 - w + a 
w - a 


H(w , a) - 


a ,/w 2 - a 2 
(w 2 - a 2 ) (a - w) 


RESULTS AND DISCUSSION 


The following four graphite/epoxy (gr/ep) laminates were used in the 
analysis: [0/145/90]^, [0], [±45] s> and [0/±45) s . Table 1 presents the 

laminate constants for the four laminates. The 0° lamina properties (ref. 5) 
were used with lamination theory to calculate the elastic constants for the 
other laminates. (Here, the 0° -direction is defined parallel to the load 
axis.) The results for the [0/±45/90] s quasi- isotropic laminate were compared 
with isotropic solutions from the literature. To show the effects of the 
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specimen boundaries, the results are presented using the following stress- 
intensity correction factors and F^: 


Kj = Sjira Fj ^11 “ $J na Fjj 

Verification 

To evaluate the improved Boundary Force Method for orthotropic materials, 
the BFM stress- intensity factor solutions were compared to results from Snyder 
and Cruse (ref. 5) for a center-crack tension specimen with a finite width of 
H/W - 3.0. In Figure 3, the curves represent the BFM calculations, while the 
symbols indicate the values taken from Snyder and Cruse. The stress -intensity 
correction factors calculated by the BFM agree, within ±3%, with the values 
from Snyder and Cruse for all laminates considered. For the [0/145/90]^ 
laminate at 2a/W - 0.8, the solution from Snyder and Cruse was 2% lower than 
the solution for an isotropic material (ref. 6), whereas the BFM solution for 
the [ 0/145/90 ] laminate was only in error by .005% compared to the accepted 
isotropic solution. Therefore, for the other laminates considered, the 3% 
difference between the present results and those of Snyder and Cruse may 
represent an error in the solution of Snyder and Cruse. 

Convergence Studies 

To demonstrate the benefit of the addition of the moment, a convergence 
study was done for a single edge crack in a quasi - isotropic plate. Figure 4 
shows the convergence of the solution for an a/W ratio of 0.6. The solution 
with both forces and moments converged to within 1.5% of the solution in ref. 
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7 with 63 DOF; the solution with only forces required 117 DOF to converge to 
the reference solution. 

Convergence studies such as that shown in Figure 4 were done for other 
laminates shown in Table 1 with an edge crack to determine if the material 
properties had any effect on the convergence rate. Although each converged to 
a different solution, the number of degrees of freedom required for 
convergence was similar for all the laminates. In general, with the moment 
included, between 80 to 100 DOF were required for convergence for 
configurations with deep cracks (a/W > 0.6). Fewer DOF are required for 
smaller a/W ratios. 


Edge Crack Solutions 

Stress- intensity correction factors were also calculated for a single 
edge crack and an inclined edge crack for a range of crack- length- to-width 
ratios and material properties. Other orthotropic solutions for these two 
problems were not available for comparison. 

Single Edge Crack . Figure 5 shows the stress- intensity correction factor 
F for a single edge crack for four graphite/epoxy laminates: [0/±45/90] s> 

[0] , [ ±4 5 ] s , and [90], The results for the quasi - isotropic laminate agreed 
within 0.005% with the isotropic results for this configuration. For very 
small a/W ratios, the solutions for the different laminates do not tend 
toward a single value, as is the case in the center crack configuration 
(Figure 3). Thus, the anisotropy of the material has more effect for small 
edge cracks than for small center-cracks. Table 2 lists correction factors 
plotted in Figure 5. 
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Inclined Edge Crack . Figures 6 and 7 show and F^, respectively, 

for an inclined edge crack for three graphite/epoxy laminates: [0/145/90] , 

[0] , and [90], In Figure 6, the mode I component does not converge to a 
single value for small a/W ratios, much as in Figure 5. The mode II 
component shown in Figure 7 varies less with anisotropy than the mode I 
component shown in Figure 6. Tables 3 and 4 list the correction factors 
plotted in Figures 6 and 7. 


CONCLUDING REMARKS 


In this paper, the Boundary Force Method (BFM) , a form of an indirect 
boundary element method, is used to analyze composite laminates with cracks. 
The BFM uses the orthotropic elasticity solution for a concentrated horizontal 
and vertical force and a moment applied at a point in a cracked, infinite 
sheet as the fundamental solution. This formulation is an improvement over 
a previous method that did not include the moment but used only the horizontal 
and vertical forces on the boundary. The necessary stress functions for this 
fundamental solution were derived using the complex variable theory of 
orthotropic elasticity. The orthotropic formulation of the BFM was verified 
by comparing solutions for a center-crack specimen subjected to uniaxial 
tension to other solutions. The BFM results agreed well with accepted 
solutions . 

Parametric studies were also done for a single edge crack, and an 
inclined edge crack, both loaded in uniaxial tension, with a variety of 
materials properties. Only a slight effect was seen due to material 
anisotropy. No other orthotropic solutions were available for comparison. 
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The additions of the moment degree of freedom to the orthotropic BFM was 
shown to greatly increase the convergence rate of the solution. Thus, fewer 
boundary elements were required to achieve the same accuracy. 

This work has resulted in a further extension of the Boundary Force 
Method in the analysis of composite laminates with cracks and notches. This 
method yields accurate solutions with minimal modeling effort, even for 
complex configurations. The accurate stress- intensity factors obtained with 
this method should be useful in predicting fracture strengths of arbitrarily 
shaped composite laminates. 
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APPENDIX A - DERIVATION OF STRESS FUNCTIONS FOR POINT LOADS 


IN CRACKED, ORTHOTROPIC SHEET 


This appendix (taken from ref. 1) presents the derivation of the stress 
function for a pair of point loads in a cracked, infinite orthotropic sheet. 
Figure 8 shows how superposition was used to determine the stress functions. 
The uncracked sheet with the point loads applied at and the stresses 

along y =» 0, |x| < a, are shown in Figure 8(b). The point loads in Figure 
8(b) are superimposed on the stresses due to the crack-face loadings in Figure 
8(c) to produce the stress-free crack face shown in Figure 8(a). 


The stress functions for a point load in an infinite, orthotropic sheet 
(used for Figure 8(b)) are derived from Lekhnitskii (ref. 2): 


^(Zj^ 



4> 2 ( z 2 ) 


B 

(z 2 - w 2 ) 


(Al) 


where 


A = C n P + C 12 Q B - C 21 P + C 22 Q 

w j “ x 0 + ^j y 0 Z j = X + ^j y (j - 1,2) 

[^(1 + % x m 2 m 2 ) + m 2 + m 2 1 

C ] i . 

2 - M 2 )(/^ 1 - ^ 1 )( M 1 - M 2 ) 
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2ni(^ l - - M 2 ) 


^ 2 K 1 + * / vx ^i^i >^2 + ^1 + ^ 
2wi(M 2 - M 1 )(a* 2 - M 2 ) (^2 ‘ 


C 


22 


M 2 [ ( M x + ^>^2 + + 

2wi(M 2 * M 1 )<M 2 * M 2 )(M 2 


^ ] 
xv 


- M X > 


(A2) 


Here v and v are the Poisson's ratios; u, . , a,, a- are 

xy yx 12 1 2 


the 


xy yx 

r 4 E X 2 ^x 1 

complex roots of the characteristic equation l/x + (— - 2i/^ )/i + — — Ol ; 

' xy xy y J 

and the barred quantities represent the complex conjugates of the underlying 
functions . 


For isotropic materials, the roots of the characteristic equation are 
- fiy - i. and the term (/^ - the denominator of equations (A2) is 

zero. When using the orthotropic equations for isotropic materials, a small 
perturbation was introduced in the values of and so that - 

i(l + «) and - i(l - e) where e - 0.0001. 

To find the stress functions for the loading shown in Figure 8(c), the 
stress functions for a point load applied at an arbitrary point on the crack 
face are used. From Savin (ref. 8) , the stress functions for this loading are 
as follows : 
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1 


^(z^ 


-p 2 p - Q 

2jt(m 1 - /i 2 > 




M 1 P + Q 1 ./a 2 - t 2 

2ff(A» 1 - m 2 ) /z 2 - a 2 Z 2 ' C 


(A3) 


where P and Q are the normal and tangential point loads applied to the 
crack face, and t is the location of the load point on the crack face (-a < 
t. < a) , as shown in Figure 9. 


By integrating equations (A3) over the crack face, the stress functions 
for the non-uniform distributed loads applied to the crack face can be written 
as follows: 


^1 ( z i ^ 


-1 


2jt(/x 1 - f. i 2 ) J z 2 - a 2 


a 
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-a 


[n 2 N(0 + T(OJ \ dg 


(A4) 


^ 2 (z 2 ) 


- p 2 ) ,/z 2 - a : 


J 


[^ X N(0 +T(0] 


J a 2 - 
z 2 - ? 


de 


In order to make the crack stress free along the crack line -a < x < a, 

the loadings N(x) and T(x) are specified to be the same as the stresses 

found in the equivalent uncracked sheet shown in Figure 8(b). That is, 

N(x) - a (x,0) and T(x) = r (x,0). Thus, from equations (3) and (Al), the 
y xy 

normal and shear stresses on the line y - 0 are 
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o (x , 0) - 2 Re (— + — ) 

y v ’ ' v x - w. x - w„' 


/^A MoB 

t (x , 0) - -2 Re ( — L + — ) 

xy v x - w 1 x - w ' 


(A5) 


Using equations (A5) and remembering that 2 Re [ f ( z ) ] = f(z) + f(z), the 
loading functions in the integrands in equations (A4) can be simplified as 
follows : 


„ 2 >.00 ♦ T(x) - 2 „ 2 Re ♦ J- 8 ^) ■ 2 Re 


(/ 1 0 - M,)A (m 2 - M]_)A (m 2 - M 2 )B 


2 *1 
x - w. 


x - x - W£ 


A R ^1 A ^ ? B 

/i.N(x) + T(x) = 2 n. Re ( — + — s ) - 2 Re ( — i + 

*1 ' ' ' ^1 x - w, x - w ' x - w. x - w„ 


(H l - /i 2 )B (n 1 - m 1 >A (a^ - 


x - w„ 


X - w. 


X - w„ 


(A6) 


Substituting these expressions into equations (A4) , integrating, and then 
adding equations (Al) , the following expressions are found for the stress 
functions for the loading shown in Figure 8(a): 
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( z i ) 


r <c 2 ■ «,)A 

[ z l - Wj f ('r ”!■ 


2(^i - ^ 2 ) / z f ' a2 


(^2 ‘ /*i )A 

a) + I f(z lt Wj, a) 

2 1 ' W 1 


(^ 2 ' /i 2 )B 

+ — f(z 1( w- ( a) + - — - 

^ • w 2 12 J Z 1 • ”l 


(A7) 


*2 <z 2> 


2 (/* L - P 2 ) ,/z 2 - a 2 


r 0^ * m 2 > b (^ - ) a 

z7 :"w. f(z 2- w 2- a) + 7“ f ( z 2’ V a) 


'2 2 


Z 2 ' W 1 


(m-j - a* 2 )b 

+ 7 — f(z 2 . w 2 , a) 

z 2 w 2 


z 2 ‘ w 2 


where f(z, w, a) - J z 2 - a 2 - Jw 7 - a 2 - z + w. 

(These stress functions are identical to those presented by Snyder and Cruse 
(ref. 5), derived by formulating the problem as a Hilbert problem.) 
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APPENDIX B - DERIVATION OF STRESS FUNCTIONS FOR MOMENT 


IN fNFINITE , UNCRACKED ORTHOTROPIC SHEET 


This appendix presents the derivation of the stress functions for a moment 
applied in an infinite uncracked orthotropic sheet. This derivation follows the 
procedure outlined in ref. 4. 


Consider the pair of forces applied as shown in Figure 10(a). By 
superimposing the stress functions for concentrated forces in the y-direction, 
the stress functions for this loading are 


^1 (Z 1 ) " C 12 P y log(z l ‘ W l ) ' C 12 P y log(z l ' W 1 ' £ > 


4 > 2 ( z 2 ) - C 22 P y log(z 2 - w 2 ) - C 22 P y log(z 2 - w 2 - e) 


(B2) 


To obtain the stress functions for a moment due to a force pair in the y- 
direction, set M = -P^e and take the limit of equations (B2) as c approaches 
zero . 


-C M 


- C 22 M 

*2< z 2> “ 


(B3) 


Consider the pair of forces applied as shown in Figure 10(b). By 
superimposing the stress functions for concentrated forces in the x-direction, 
the stress functions for this loading are 
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^ 1 (z 1 ) - C 11 P x log(z 1 - w x ) - C 11 P x log(z 1 - w 1 + 


<t> 2 ( z 2 ) 


C 21 P x log(z : 


- w 2 ) - 


C 21 P x log(z 2 


- w 2 + /i 2 e y ) 


(B4) 


To obtain the stress functions for a double force due to a moment in the x- 
direction, set M = *P X « and take the limit of equations (B4) as e approaches 
zero . 


<^ 1 (z 1 ) 


Z 1 - W 1 


^ 2 <z 2 ) 


p 2 C 21 M 
z„ - w„ 


(B5) 


By superimposing these stress functions, the solution for a concentrated 
moment due to the two force pairs (M - -2Pe) can be written as 


. . , ( "i c n - c i?>" 

W ' 2 (z x - *j) 


' -(^, c n - c )M 

W - " 2 }, - 


i 2 ( z 2 ) 


( / *2 C 21 ' C 22 )M 
2(z 2 - w 2 ) 


*2 (z 2* = 2 (z 


- (^2 C 21 ~ C 22 )M 

w 2 ) 2 


(B6) 


(B7) 
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Table 1 - Laminate Constants 


Laminate 

E , GPa 

X 

E , GPa 

y 

V 

xy 

G , GPa 
xy 

[ 0/±45/90 ] s 

60.04 

60.04 

0.259 

23.85 

[0] 

11.72 

144.8 

0.017 

9.65 

[ ± 45] s 

31.18 

31.18 

0.615 

38.03 

[±30] s 

17.64 

66.99 

0.262 

30.94 

[90] 

144.8 

11.72 

0.210 

9.65 
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Table 2 - Stress - Intensity Correction Factors for Edge Crack 


a/w 


F 




[0/±45/90] s 

[0] 

[±45] s 

[90] 

.1 

1.1894 

1.1589 

1.2861 

1.1889 

.2 

1.3656 

1.3576 

1.4537 

1.3600 

.3 

1.6567 

1.6923 

1.7623 

1.6397 

.4 

2.1062 

2.1884 

2.2088 

2.0594 

.5 

2.8015 

2.9034 

3.0347 

2.9772 

.6 

4.0021 

4.2344 

4.2323 

3.7939 
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Table 3 - Stress - Intensity Correction Factors for Inclined Edge Crack 


a/w 


F I 



[ 0/145/90 ] g 

[0] 

[90] 

.1 

0.7113 

0.7495 

0.7164 

.2 

0.8130 

0.8180 

0.8043 

.3 

0.8713 

0.8987 

0.8759 

.4 

1.0091 

1.0407 

1.0113 

.5 

1.1952 

1.2095 

1.1969 

.6 

1.4442 

1.4792 

1.4462 
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Table 4 - Stress - Intensity Correction Factors for Inclined Edge Crack 


a/w 


F II 



[0/±45/90] s 

[0] 

[90] 

.1 

0.3739 

0.3737 

0.3760 

.2 

0.3957 

0.4029 

0.4001 

.3 

0.4452 

0.4451 

0.4472 

.4 

0.5024 

0.5026 

0.5047 

.5 

0.5795 

0.5804 

0.5815 

.6 

0.6804 

0.6855 

0.6820 
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Superposition for formulati 
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Figure 4. Comparison of convergence of solution with forces only and with 
forces and moments for single edge crack. H/W - 2.0, a/W - 0.6 




ge crack. 








Figure 9. Point loading on crack face. 
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